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ABSTRACT 

Nitrogen is the fifth most abundant element in the Universe, yet the gas-phase chemistry of N-bearing species remains poorly understood. Nitrogen 
hydrides are key molecules of nitrogen chemistry. Their abundance ratios place strong constraints on the production pathways and reaction rates 
of nitrogen-bearing molecules. We observed the class protostar IRAS 16293-2422 with the heterodyne instrument HIFI, covering most of the 
frequency range from 0.48 to 1.78 THz at high spectral resolution. The hyperfine structure of the amidogen radical 0-NH2 is resolved and seen in 
absorption against the continuum of the protostar. Several transitions of ammonia from 1.2 to 1.8 THz are also seen in absorption. These lines trace 
the low-density envelope of the protostar. Column densities and abundances are estimated for each hydride. We find that NH:NH,:NH 3 a;5:l:300. 
Dark clouds chemical models predict steady-state abundances of NH 2 and NH 3 in reasonable agreement with the present observations, whilst that 
of NH is underpredicted by more than one order of magnitude, even using updated kinetic rates. Additional modelling of the nitrogen gas-phase 
chemistry in dark-cloud conditions is necessary before having recourse to heterogen processes. 

Key words. ISM: abundances, Astrochemistry, ISM individual objects: IRAS 16293-2422 
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1. Introduction 

Nitrogen is the fifth most abundant element in the Universe 
and is a fundamental component of molecules associated with 
life. Nitrogen-bearing molecules are routinely observed to- 
wards a wide vari ety of environments, from the diffuse in- 
terstellar medium ( Liszt & Lucasl |2001|) to pre- stellar cores 
dBergin & Tafallall2007h and protoplanetary disks dDutrev et al.l 
Il997l) . Complex N-b earing molecules are also obser ved towards 
star-forming regions dHerbst & van Dishoec]dl2009l) . The chem- 
ical network of nitrogen is apparently simple in that a small set 
of reactions is involved (Pin eau des Forets et al.lll990 l hereafter 
PdF90). The chemistry of nitrogen has been modelled in vari- 
ous environments with moderate suc c ess, including ty pical dark 
cloud conditions (iMillar et al.l 1 1 99 lj: | L e Bourlo till99ll). shocks 
(PdF90), pre-stellar co res dFlower et al.ll2006tlMaret et al.f2 006 ; 
iHilv-Blant et al.l 2010|). and photo-dissociation regions (PDR) 
Sternberg & Dalg arnol[T995l) . One major unknown is the total 
abundance, in dense and shielded environments, of gas-phase ni- 
trogen, the reservoir of which consist of N and/or N 2 . Because 
they are not directly observable, estimates of their abundances 
rely on observations of other N-bearing compounds and chemi- 
cal modelling. Observational constraints of the dominant chemi- 
cal pathways of the nitrogen chemistry and their kinetic rates are 
thus crucial. 



In this respect, nitrogen hydrides are of utmost importance 
since they are among the first neutral N-bearing molecules 
formed in an initially atomic gas dominated by hydrogen and 
helium. Ammonia was among the first interstellar molecules de- 
tected in emission towards the Galactic centre ( Cheun sTet al.l 
1 19681) . The li ghtest radical, i midog en NH, was observed in 
absorption by iMever & Rothl dl99lf) along the diffuse line of 
sight towards £Per. Gas-phase models were found to underes- 
timate the abundance of NH, and dust gr ains were then pro- 
posed to solve part of the discrepancy ( Mey er & Rothl 1 1 99 1 : 



IWagenblast et al.l Il993l). Amidogen (NF L) was observed in 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led principal Investigator consortia and with impor- 
tant participation from NASA. 



absorption by Ivan Dishoeck et al. I (fl993l) fr om dense gas in 
Sgr B 2. This source was also targeted by iGoicoechea et al.l 
(2004) with ISO. The ratios of the three hydrides were found 
to be NH:NH 2 :NH 3 s! 1:10:100, incomp atible with t h e dar k 
cloud value NH 3 /NH 2 < 3 predicted by IMillar et all dl991l) . 
Unfortunately, the modelling of the chemistry in Sgr B2 is dif- 
ficult due to the complexity of the source, which probably in- 
corporates shock dynamics. The NH:NH 2 :NH 3 ratios measured 
in Sgr B2 may thus not be representative of cold dark clouds, 
and the chemistry of nitrogen hydrides in these environments re- 
mains largely unexplored. 

In this paper, we present HIFI observations of the submil- 
limetre lines of NH 2 and NH 3 in absorption against the contin- 
uum of the class protostar IRAS 16293-2422. Section 2 sum- 
marizes the observation strategy and data reduction. In Sect. 3 
we derive the column densities of NH 2 and NH 3 . Abundances of 
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the three lightest nitrogen hydrides are estimated in Sect. 4, and 
compared to steady-state models using updated reaction rates. 

2. Observations and data reduction 

The solar-mass protostar IRAS 16293-2422 was observed with 
the HIFI instrument onboard the Herschel Space Observatory, 
as part of the HIFI g uaranteed time key program CHESS 
dCeccarelli et al.l [2010). Full spectral coverage of bands la 
(480 - 560GHz), 3b (858 - 961 GHz), 4a (949 - 1061 GHz), 5a 
(1.12 - 1.24 THz), and 7a (1.70 - 1.79 THz) were performed on 
2010 Mar 1, 3, and 19, using the Spectral Scan DSB mode with 
optimization of the continuum. The Wide Band Spectrometre 
(WBS) was used as a backend, providing us with a spectral reso- 
lution of 1 . 1 MHz over an instantaneous bandwidth of 4x 1 GHz. 
The targeted coordinates were £#2000 = 16 /! 32'" 22'75, £2000 = - 
24° 28' 34.2". The two reference positions were situated approx- 
imately 3' east and west of the source. The beam size is well 
approximated by HPBW = 21.5"/vthz- For the analysis, inten- 
sities were then brought to a main-beam temperature scale using 
Feft = 0.96 and B eS = 0.70. 

The d ata were processed using the standard HIFI pipeline 
HIPE 2.8 (Ott 2010) up to frequency and amplitude calibrations 
(level 2). For the SIS bands 1 to 5, a single local-oscillator- 
tuning spectrum consists of 4 sub-bands of 1 GHz for each polar- 
ization. The 1 GHz chunks for bands 1 to 5 are then exported as 
FITS files in the CLASS90/GILDAS formafl dHilv-Blant et al l 
120051) for subsequent data reduction and analysis. Despuring 
and residual bandpass effect subtraction were performed in 
CLASS90 using generic spectral-survey tools developed in our 
group. Sideband d econvolution is computed with the minimiza- 
tion algorithm of IComito & Schilka (l2002h implemented into 
CLASS90. Line identification used the Weeds CLASS90-add- 
on developed bv lMaret et al.l d2010l) . which provides an efficient 
interface to the public CDMS and JPL spectroscopic databases 
(Mul ler et alJ200^lPickettetd]IT998h . 

3. Results 

Figure Q] shows the detection of the hyperfine structure (HFS) of 
the N - - 1, J — I - |, and / = \ — A transitions of amidogen 
in its ortho form, with their strongest components at the rest fre- 
quencies 952 5783 54 GHz and 959.511716 GHz, respectively 
dMtiller etal.lll999l) . The HFS is almost entirely resolved with 
an intensity ratio that clearly deviates from optically thin LTE 
excitation. The 3cr noise levels are indicated in each panel. The 
para-NH 2 line was not detected and a 5cr upper limit to the main 
HFS component at 947.725 GHz is J t dv < 0.2 km s" 1 . Figure[2] 
shows several transitions of ammonia seen in absorption, from 
1168.4 to 1763.8 GHz. The fundamental rotational transition 
at 572.6 GHz was also detected in emission (see Fig. IA.31 > but 
because it contains both emission and absorption, it is not dis- 
cussed in this paper. 

All column densities are derived assuming a single excita- 
tion temperature r ex for each molecule. The opacity of each HFS 
component and the excitation temperature are determined by si- 
multaneously fitting all HFS components, constraining the opac- 
ities to scale with A u ig u . The fit is performed in CLASS90 by 
applying the HFS method to the continuum-subtracted spectra 
T m b = [J v (T ex ) - J v (T C mb) - T c ] (1 - e~ T "0- The total column 
density reads AT tot = (8ttv 3 /c 3 ) Q/(A ulgu ) x / t uI dvxe E <' kT " /[I - 

1 http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 1. Absorption spectra of o-NH 2 at 952 GHz and 959 GHz, 
with the HFS fit overlaid (line). The relative intensities in the 
optically thin LTE limit are indicated at the bottom. 

e -hv/kT cx ^ w here Q is the partition function at T ex . In deriv- 
ing the excitation temperature, we assumed equal filling fac- 
tors for the absorbing gas and the background continuum radi- 
ation. The SSB continuum intensity Tc, needed to derive r ex , 
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Fig. 2. NH 3 absorption lines from 1.168 to 1.764 THz. LTE pre- 
dictions are shown in red (r ex =10 K, /V(NH 3 ) = 3.5xl0 15 cm" 2 ) 
and blue (r ex = 8 K, /V(NH 3 ) = 2xl0 16 cnr 2 ). 
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Table 1. Column densities of nitrogen hydrides towards IRAS 16293-2422. 



Species 


Transition 


C omponent fl 


Frequency 


HPBW 




Tf 






FWHM" 


Mot 








GHz 


arcsec 


K 


K 




K 


km s -1 


lO'W 2 


Nff 


Oi -» 1 


2 ' 2 2 ' 2 


946.476 


23 


0.8 


-3.8 + 0.6 


2.8 ± 0.7 


9.5 


0.60 


2.20 ± 0.80 


o-NH 2 


°ooi -» ln| 


3 5 5 7 

2 ' 2 2 ' 2 


952.578 


23 


0.9 


-9.0 + 0.5 


12.8 + 0.7 


8.5 


0.60 


0.40 ± 0.06 




°ooi -> l u\ 


3 5 3 5 
2' 2 2' 2 


959.512 


22 


0.9 


-2.5 ±0.2 


4.9 + 0.7 


9.5 


0.60 


0.59 ±0.12 


p-NH 3 


1-2 


(!,-)-(!,-) 


1168.453 


18 


1.2 




300 - 70 


8-10 


0.50 


200 - 35 


o-NH, 


1-2 


(0, +) - (0, +) 


1214.853 


18 


1.3 




470 - 130 


8-10 


0.50 


200 - 35 


p-NH, 


1-2 


(l,+)-(l,+) 


1215.246 


18 


1.3 




330 - 80 


8-10 


0.50 


200 - 35 


p-NH 3 


2-3 


(2,-) -(2,-) 


1763.823 


12 


2.2 




2.0 - 1.4 


8-10 


0.50 


200 - 35 



Notes. {a) For NH, the quantum number s for the rotational transition Nj are Fx = I H + J and F = I N + ¥i i Klaus et al. 1997). For the N KaKcJ rotational transition of NH 2 , the quantum 
numbers are Fi = I N + J and F = I H + Fi i Miiller et al. 199§. In the case of ammonia, quantum numbers are given separately for J = N + S and (K, e), where 6 is the symmetry index (see 
Maret et al. 2009). For these lines, the frequency given is that of the brightest HF component. Single-sideband continuum in a T m b scale. * c) 77 = r B /[./y(7' ex ) - J v (Tcmb) - TcmA from 
the HFS fit. In the case of ammonia, see Sec t. 3. ^ A conservative uncertainty of 0.25 MHz (0.08 km s ) imposed by the HIPE 2.8 pipeline was retained. '''' The integrated line opacity is 
calculated as /rdv = 1.06 F WHM X t,,,. w IBacmann et aljt20T(l) . 



is estimated as half the median of each 1 GHz chunk computed 
in line-free spectral windows prior to deconvolution, assuming 
equal gains in the two sidebands (see details in Appendix [A}. 
It is found to increase with frequency and, for rest frequen- 
cies vthz — 0.492 to 1.242 THz, Tc is well approximated by 
7c(Tant, [K]) = 1.10 vthz - 0.42 (see Fig. lAlll . At higher fre- 
quencies, Tc is estimated from the deconvolution tool in the 
HIPE software, assuming equal sideband gains. 

For ammonia, the hyperfine components in each lines are not 
resolved out by our observations, so the column density and ex- 
citation temperature were determined by simultaneously fitting 
the 4 lines shown in Fig. [2] The 1763.8 GHz line has a well con- 
strained centre-line opacity, yet the excitation temperature and 
column density remain degenerate. We therefore varied the exci- 
tation temperature between 5 to 15 K, and adjusted accordingly 
the column density until the observed line profiles were repro- 
duced. The line width was fixed to 0.5 kms 4 . For T ex > 10 K, 
the absorption lines at 1.11684, 1.2148, and 1.2152 THz become 
significantly weaker than observed. In contrast, for T ex < 8 K, 
the column density needed to reproduce the 1.7 THz line be- 
comes so large that the 1.2 THz lines broaden significantly. 
The full set of constraints points towards T ex as 8 - 10 K and 
Mot ~ 2.0xl0 16 - 3.5xl0 15 cirr 2 . Figure shows the results of 
the modelled lines corresponding to T ex = 8 and 10 K. The non- 
detection of the 1.7635 and 1.7636 THz absorption lines (not 
shown in Fig. |2j is consistent with these models. 

The results are summarized in Table Q] The excitation tem- 
perature of the 952 and 959 GHz transitions of o-NH 2 are sig- 
nificantly higher than 2.73 K. The different values for each HFS 
most likely results from the LTE assumption not being entirely 
valid. However, T ex <k hv/k, so that the corresponding uncer- 
tainties have negligible consequences on the column density. 
The thermalized ortho:para ratio is expected to be large at low 
temperatures (« e 3039 ^ T , see Fig. IA.2I >. Our two determinations 
of the o-NH 2 column density thus give an average for the total 
NH 2 column density. The resulting NH 2 and NH 3 column den- 
sities are 0.44 + 0.07xl0 14 cirr 2 and * 35 - 200xl0 14 cirr 2 , 
respectively. The column density ratios are thus NH/NH 2 = 
5.0 ± 1.2 and NH 3 /NH 2 = 90 - 500, or NH:NH 2 :NH 3 «5: 1:300. 



4. Discussion 

From Table [1] the low excitation temperatures indicate that all 
lines are sub-thermally excited and therefore most likely trace 
regions with densities much lower than their critical densities, 
which are of the order of 10 7 cirr 3 . The physical source model 



of Crimie r et al.l d2010h predicts densities lower than « 10 6 cirr 3 
for radii larger than 2400 AU or 17" at a distance of 120 pc, 
comparable to the HPBW of the present observations, thus sup- 
porting the assumption of equal filling-factors for the absorbing 
gas and the continuum emission. At these radii, the modelled 
gas temperature is lower than 20 K. The NH and NH 2 linewidths 
are thus dominated by non-thermal broadening suggesting that 
turbulence has not been dissipated in the absorbing gas. 

To place constraints on the average abundances in the fore- 
ground absorbing material, we need to measure the total H 
column density, which is no t directly observabl e. The one- 
dimensional density profile of ICrimier et alJ d2010h is extended 
to lower densities as n Hl (r) = 3xl0 8 (r/85AU)" 18 + 10 3 cm 4 , 
to allow for a low-density envelope. The column density profile, 
convolved by the HIFI beam, is dominated by the lines-of-sight 
close to the centre. Considering only the gas in the regions with 
«h 2 < 10 6 cm 4 , the column density is 7V(H 2 ) ~ 8xl0 22 cirr 2 
or 80 magnitudes of visual extinction (assuming standard dust 
properties and that all H is molecular). Were the hydrides ab- 
sorption to occur at densities lower than 10 4 or 10 5 cm 4 , the 
total column density would decrease by a factor 3. Therefore, 
in the following we estimate the abundances assuming 7V(H2) = 
5.5 + 2.5xl0 22 cm 4 . Results are summarized in Table|2] 

The column density ratios observed in th e cold envelope of 
IRAS 16293-2422 differ from those derived bv lGoicoechea et aD 
(2004]), who found NH:NH 2 :NH 3 = 1:10: 100 towards Sgr B2. 
These authors noted that these ratios are not consistent with typ- 
ical dark cloud conditions but can be explained by shock chem- 
istry. Though the NH 3 /NH 2 ratios are similar in both the cold en- 
velope and the Sgr B2 region, the NH/NH 2 ratios are drastically 
different, suggesting that different chemistry is at work in the 
two sources. To study this, we have employed the A strochem 
gas-phase chemical code of Maret & Bergin (2010) combined 
with modified versions of the osu . 09 . 2008 chemical networifl 
Typical physical conditions for a cold molecular cloud (gas tem- 
perature T — 10 K, «h = 2 x 10 4 cm 4 , A v = 10 mag) and 
a cosmic -ray ionization rate ( - 1.3 x 10 47 s 4 were used. 
Higher extinctions would not modify the predicted steady-state 
abundances. In a similar way, higher densities would shorten 
the time to reach a steady-state but would not alter the cor- 
responding abundances. T he initial abundances are taken from 
IWakelam & Herbstl (120081) for which the C/O gas-phase abun- 
dance ratio is 0.41. The resulting time-dependent fractional 
abundances of the nitrogen hydrides are shown on Fig. IB. II 
The steady-state is reached after a few 10 6 yrs, for which the 



2 http : //www . physics . ohio-state . edu/~ eric/research . html 
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Table 2. Predicted steady-state fractional abundances of several 
nitrogen-bearing species in three different models. 



Species 


I 




II 




III 




Observations 


N,H +D 


100:0 


90:10 


90:10 




NH+'' 


85:2: 


13 


85:2: 


13 


95:2: 


3 




N 2 


1.0(- 


■5) 


1.0(- 


■5) 


1.0(- 


-5) 




N 


2.4(- 


-7) 


2.5(- 


-7) 


2.9(- 


-7) 




NH 


3.5(- 


-11) 


9.5(- 


■11) 


8.8(- 


-11) 


2.0 ± 1.0(-9) 


NH, 


1.9(- 


■10) 


1.9(- 


■10) 


1.3(- 


-10) 


4.0±2.0(-10) 


NH~ 


3.7(- 


-8) 


3.7(- 


-8) 


8.3(- 


-8) 


0.3^! -2tl(-8) 


NH:NH,:NH, 


0.2:1 


:190 


0.5:1 


:190 


0.7:1 


:640 


5:1:300 



Notes. Numbers in parenthesis are powers of ten. Branching ratios (BR) for the disso- 
ciative recombination (DR) of NH+ are N + H 2 (4%), NH + H (39%) and N + 2H (57%) 
IThoma s et al. 2005). <o) BR for the DR of N,H + into N, + H and NH + N, respectively 
i Molek et al. 2007tlAdams et alJ20091) . <6) BR for the DR of NH+ into NH 3 + H, NH 2 + H 2 , 
and NH 2 + 2 H, respectively jOiekull etTal . 2004). See additional details in Table lB.ll 



dominant chemical paths (dissociative recombinations, neutral- 
neutral reactions at low temperature) shall also be explored in 
the process. 
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NH:NH 2 :NH 3 abundance ratios are 0.2:1:190 (see TableEJ. We 
note that the NH:NH 2 :NH 3 ratios are roughly constant for times 
> 10 5 yrs. If NH 3 /NH 2 is consistent with the observations to 
within a factor of 2, the steady-state NH/NH 2 ratio is too small 
by more than one order of magnitude because the abundance of 
NH is underpredicted by our model. 

NH 2 and NH 3 are princip ally formed by th e dissociative 
recombination (DR) of NH+ dLe Bourlotl 1199 ll) . The forma- 
tion of NH is dominated by NH 2 (0,OH)NH, unless the NH 
channel of the DR of N 2 H + has a non-zero branching ra- 
tio (BR). There are still disparate results for the BR of these 
PR reactions, and the most recen t literature (see Table |2]and 
Florescu-Mitc hell & Mitchell l2006h suggest uncertainties of at 
least 10% in the BR. We conducted three model calculations, 
where we varied the BR by 10%, to explore the effects on 
NH:NH 2 :NH 3 . The results are summarized in Table|2] The abun- 
dances vary at most by factors of a few. Noticeable is the increase 
in «(NH) associated with the opening of the NH channel of the 
DR of N 2 H + . In any case, NH remains underabundant by al- 
most two orders of magnitude. This deficit is reminiscent of the 
early results obtained in the diffuse medium that motivated the 
recourse to surface reactions. We recall that our models do not 
include freeze-out onto dust grains to enhance the role of gas- 
phase reactions. In the case of NH formation, however, the role 
of surface reactions is poorly constrained. 

5. Conclusions 

We have presented absorption spectra of the hyperfine struc- 
ture of NH 2 and several transitions of NH 3 . These lines em- 
anate from the low density envelope of the protostar, at den- 
sities lower than typically 10 6 crrT 3 , representative of typical 
dark cloud conditions. We have determined the column den- 
sities, as well as fractional abundances. We have found that 
NH:NH 2 :NH 3 ^5:1:300. On the basis of on an updated chemi- 
cal network, we have computed the steady-state abundances of 
amidogen and ammonia, which agree well with observed values. 
Imidogen is undepredicted by more than one order of magni- 
tude, and NH/NH 2 < 1. At this point, more observations (of e.g. 
15 N isotopologues) and modelling are clearly needed. Freeze- 
out and surface reactions, not included in our calculations, are 
potentially important. Before having recourse to dust surface 
processes, however, gas-phase chemistry has to be explored in 
far more d etails in studying e.g., t he influence of the gas-phase 
C/O ratio dHilv-Blant et all 12010) or the effects of the H 2 o:p 
ratio. The consequences of the uncertainties in the rates of the 
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Fig.A.l. Double sideband continuum (in antenna temperature 
scale) level towards IRAS 16293-2422, from 0.5 to 1.2 THz. The 
straight line is the result of a linear fit. The observed HIFI bands 
are highlighted. 



Appendix A: Determination of the column density 



Figure I A. 1 1 displays the double-sideband continuum intensity 
measured from 0.5 to 1.2 THz with the double-sideband re- 
ceivers of the HIFI instrument. The single-sideband continuum 
is then estimated assuming equal image and signal gains. The 
SSB continuum intensity was estimated following two indepen- 
dent methods. First, the median of the intensity in each 1 GHz 
sub-band was computed. The resulting double-sideband (DSB) 
continuum level Tq, sampled every 0.25 GHz from 0.492 THz 
to 1 .242 THz, was found to increase linearly with the frequency 
(see Fig. IA.U . The second method consisted in deconvolving 
the despured spectra prior to baseline subtraction, the result of 
which is the SSB continuum level. The DSB deconvolution was 
applied to the subset of spectra covering the spectral ranges of 
the lines considered in this paper. The DSB continuum level is 
twice the SSB to better than 10% in these intervals. From the first 
method, the antenna temperature scale SSB continuum level is 
wellfittedby a 1st order polynomial as Tc[K] = 1.10 vthz _ 0.42, 
for rest frequencies vjh z ranging from 0.492 to 1.242 THz (see 
Fig. IA.U . The increase in Tc with frequency might be caused 
by the convolution of the dust temperature profile with the tele- 
scope beam. It may also trace the increase in the dust emissivity 
with frequency. At the higher frequencies of the ammonia lines 
in Band 7, the continuum level was estimated using the second 
method, in HIPE. 

Figure lA~2l shows the o:p ratio for NH 2 assuming equilibrium 
at a single temperature for all levels. 

Appendix B: Chemical modelling 

Nitrogen chemistry starts with the formation of NO and CN by 
means of the re actions N(OH,NO)H a nd N(CH, H)CN, which 
then lead to N 2 (iHilv-Blant et all 1201 0l) . Once N 2 is formed, it 
reacts with He + to form N + , which, by successive hydrogen ab- 

and NH4. The dis- 



stractions, leads quickly to NH , NH 2 , NH, , 



sociative recombination (DR) of NH| is the dominant forma- 
tion route for NH 2 and NH 3 . The formation of NH is dominated 
by NH 2 (0,H)NH, unless the NH channel of the DR of N 2 H + 
has a non-zero branching ratio (BR). Although DR branching 
ratios have been (re)measured recently for N 2 H + , NH 2 , and 
NH4, there are still disparate results and significant uncertain- 
ties. According to the most recen t litterature (see Table |2]and 
Florescu-Mitchell & Mitchell 2006), the BR of DR reactions are 
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Fig. A.2. Thermalized ortho:para ratio for NH 2 . 
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Fig. A.3. Emission line of the fundamental transition of NH 3 
572.6 GHz detected by HIFI towards IRAS 16293-2422. 



at 



uncertain by at least 10%. We conducted three model calcula- 
tions, where we varied the BR by 10%, to explore the effects on 
NH:NH 2 :NH 3 . The results are summarized in Table|2] The abun- 
dances vary at most by factors of a few. We note the increase in 
n(NH), which is associated with the opening of the NH channel 
of the DR of N 2 H + . In any case, NH remains underabundant by 
almost two orders of magnitude. This deficit is reminiscent of 
the early results obtained in the diffuse medium that motivated 
the recourse to surface reactions. We recall that our models do 
not include freeze-out onto dust grains to ensure that the role of 
gas-phase reactions is enhanced. In the case of NH formation, 
however, the role of surface reactions is largely ill-constrained. 

To model the steady-state abundances of the nitrogen hy- 
drides, the osu. 09 . 2088 network and rates were used. The 
osu.09.2008 network contains 13 elements, 449 species, and 
4457 gas-phase reactions. We note that this version of the OSU 
database does not contain molecular anions or any depletion of 
gas-phase species. The network has been updated from a revi- 
sion of the branching ratios and rate coefficients for the dissocia- 
tive recombination (DR) of the nitrogen bearing cations N 2 H + , 
NH 2 , and NH4: 
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Fig. B.l. Predicted abundances for several N-bearing molecules 
as a function of time for model 1 ("standard", see text), for a 
density «h = 2xl0 4 crrT 3 , a gas temperature 7\i n =10 K, and a 
total extinction Ay=10 mag. Abundances are given with respect 
to H nuclei. 



Table B.l. Dissociative recombination branching ratios em- 
ployed in the present work. 



Ion Total rate" Products Model 

y /? 1 2 3 

N 2 H +fc 1.00(-7) I 03 N, +H 100% 90% 90% 

NH + H 0% 10% 10% 

NHf 3.00(-7) -0.5 N + H + H 57% 

N + H 2 4% 

NH + H 39% 

NHj 3.10(-7) -0.5 NH + H + H 50% 

NH, + H 50% 

NH+ rf 9.40(-7) -0.6 NHj + H 85% 85% 95% 

NH^ + H, 2% 2% 2% 

NH^ + H + H 13% 13% 3% 



Notes. M Reaction rates are written as k = ■yf T/iOOf in cm'r 1 . " J) IMolek et all 
i 2007)«'> Thomas et al. 1 2005) <<" Ojekull et al. 1 2004). 



- The DR of N 9 H + has been determ ined experimentally us- 
ing both flowing afterglow (F A, see lAd ams et al. 1 199 lb and 
storage ring (SR) techniques dGeppert et alj20 04). leading to 
controversial results concerning the branching ratios (BR) of 
the two channels N 2 + H and NH + N. Thus, in contrast to the 
FA results that established the major product as N 2 + H with 
a BR « 100%, iGeppert et all (120041) found this channel to 
account for only 36% of th e total reaction. The most recent 
FA an d SR measurements dMolek et alj|2007t lAdams et all 
2009), however, have confirmed the earlier FA results that 
the DR of N 2 H + should lead predominantly to N 2 + H with 
a BR ss 90 - 100%. This result is also supported by the ab 
initio calculations of iTalbil (l2009h . For the total rate coeffi- 
cient, we adopted the (temperature dependent) expression of 
IGeppert et all d2004 . as in the osu. 09. 2008 network. 

- For the DR of NH^ , the latest SR measurements are those of 
iThomas et al] d2005h . who obtained the BR for N + H 2 (4%), 
NH+H (39%), and N+2 H (57%). For th e rate coefficient we 
adopted the expression recommended by Mi tchelll dl990h . as 
in the osu. 09 . 2008 network. 

- Finally, the SR measurements of lOiekull et al] d2004b 
demonstrated that the DR of NH^ is dominated by the 
product channels NH 3 + H (85%), NH 2 + 2H (13%), and 
NH 2 + H 2 (2%). For the total rate, we adopted the expression 
of Oi ekull et al.l d2004b . which differs slightlty from the one 
recommended in the osu . 09 . 2008 network. 

We note that for NH3 there is to our knowledge neither measure- 
ments nor calculations available. We therefore adopted the rate 
and branching ratios recommended in the osu. 09. 2008 net- 
work, corresponding to the two channels NH + 2H (50%) and 
NH 2 + H (50%). 

Following the above update of DR rates and branching ra- 
tios, we have adopted three different chemical models by vary- 
ing the DR branching ratios of N 2 H + and NH4 within 10%, the 
typical experimental uncertainty. The DR of these two ions was 
indeed found to be the dominant formation routes of NH and 
NH 2 , respectively. The employed branching ratios are listed in 
TablelO 



